We sought to determine the scale of scintillation in the in- 
terstellar plasma of PSR J0437-4715. We used the Very Long 
Baseline Array to obtain scintillation amplitude and phase data, 
from dynamic spectra at 327 MHz. We observe two scales of 
scintillation of pulsar PSR J0437-47 15, differing by more than 
an order of magnitude in scintillation bandwidth. The wider- 
bandwidth scale of scintillation that we observe indicates less 
scattering for this pulsar than for other nearby pulsars, other 
than PSR B0950+08. 
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Abstract. 



1. Introduction 

In this paper we report measurement of the interstellar scatter- 
ing by plasma along the line of sight to pulsar J0437-4715. 
This is one of the closest pulsars. It has distance R = 150 pc, 
and its transverse velocity is 100 km/sec (van Straten 200T)- It 
is quite strong over a wide range of observing frequencies. Its 
proximity, and its intensity, make it an important probe of the 
local interstellar plasma. 

We describe observations of the scattering parameters of 
J0437-4715, made at 327 MHz. We find two scales of scin- 
tillation, corresponding to two scales of structure in the pul- 
sar's dynamic spectrum. Multiple scales are not uncommon for 
nearby pulsars (see, for example, Stinebring et al. 2001 Hill 
et al. 2003 1, although scales separated by more than an order 
of magnitude are unusual. Comparison with other pulsars of 
similar distance and dispersion measure shows that this pulsar 
has relatively less scattering than expected on the basis of a 
model for homogeneous turbulence. In this paper, we discuss 
the observations and comparison with other pulsars. We dis- 
cuss comparison with previous observations of scintillation of 
J0437-4715, and interpretation of the measurements, in an ac- 
companying paper (Smirnova et al. 2006 1, hereafter Paper 2). 

2. Observations and Analysis 

2.1. Data Collection and Calibration 

We observed PSR J0437-4715 at two epochs using the Very 
Long Baseline Array (VLBA) operated by the National Radio 
Astronomy Observatory (NRAO). On 5 November 1996, we 
observed the pulsar over four 33-minute scans, using 25-meter 
telescopes at Fort Davis, Mauna Kea, Pie Town, and St. Croix, 
with a 32-MHz bandwidth centered at 332 MHz, and with a 

* Observations made with the US VLBA. The National Radio 
Astronomy Observatory is a facility of the National Science 
Foundation operated under cooperative agreement by Associated 
Universities, Inc. 



250-kHz spectral resolution. Between the scans, we observed 
extragalactic sources for calibration purposes. On 8 April 1999, 
we observed the target for two 44-minute scans and one 33- 
minute scan, using the telescopes at Fort Davis, Kitt Peak, Los 
Alamos, Mauna Kea, Owens Valley, Pie Town, and St. Croix, 
with a bandwidth of 32 MHz centered at 328 MHz, and with a 
spectral resolution of 125 kHz. Mutual visibility is short, par- 
ticularly on long baselines, because the source lies at low decli- 
nation; this limited the useful time span of the observations. To 
increase the signal-to-noise ratio, the data for this observation 
were correlated using a gate which covered 12% of the pulse 
period and which isolated the central peak of the pulse profile. 
Scans of a calibrator source preceded and succeeded the obser- 
vations of the pulsar. For both observations, we aggregated the 
32-MHz bandwidth using four 8-MHz intermediate-frequency 
(IF) bands. 

We performed the first half of the calibration using NRAO's 
Astronomical Image Processing System (AIPS). The data from 
both observations were marred by radio-frequency interference 
(RFI). The second observation in particular suffered badly from 
RFI, forcing us to eliminate the Kitt Peak and Owens Valley 
antennas from the outset. The interference also tainted some 
measurements of system temperature. The system temperature 
should not change abruptly, so if an isolated reading clearly de- 
parted from an otherwise smooth trend, we removed that value. 
In two IFs, however, rampant corruption due to RFI obliter- 
ated evidence of a smooth evolution in system temperature. For 
these two IFs, we copied the system temperatures from an ad- 
jacent IF. 

We examined the remaining visibilities and flagged chan- 
nels and time ranges to eliminate data corrupted by very strong, 
narrow-band signals typical of man-made sources. We applied 
the usual corrections for the effects of antenna parallactic an- 
gles, digital sampler bias, system temperature, and gain calibra- 
tion. We then performed a manual phase calibration by running 
a fringe fit on a thirty-second interval of calibrator data to deter- 
mine the electronic delay differences between IFs, and we re- 
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moved the resulting phase differences from the pulsar and cal- 
ibrator data. The final procedure done using AIPS was a fringe 
fit on each source to remove phase slopes due to the residual 
fringe rate and delay. 

To perform a bandpass calibration, we averaged, for each 
baseline, the complex visibilities from a ten-minute interval 
of calibrator data. Deviations from a fiat amplitude and phase 
spectrum reflected the non-ideal frequency response of the 
baseband converters. The gain of the filters drops off at the ends 
of each IF, significantly reducing the signal-to-noise ratio in 
the first and last channel of each IF; consequently, we omitted 
these channels in our subsequent analysis. We then smoothed 
the data using boxcar averaging with a one-minute window to 
further increase the signal-to-noise ratio. We then exported the 
data from AIPS and continued calibration using custom soft- 
ware. 

At our observing frequency, propagation through the iono- 
sphere significantly affects the visibility phases in both time 
and frequency (Thompson et al. 1986 1. Data from both obser- 
vations exhibited broadband variations in phase of order half a 
radian over a time scale of about fifteen minutes, as would be 
expected for ionospheric path changes. We approximated and 
removed these temporal variations using an eighth-order poly- 
nomial fit to amplitude- weighted phase averages calculated for 
each time point. 

The ionosphere also introduces a nonlinear variation of 
phase with frequency (Thompson et al. 1986 >. The fringe 
fit removes only a linear dependence of phase on frequency. 
Usually, the bandwidth is small enough compared to the ob- 
serving frequency that the higher-order terms are not impor- 
tant. However, in our case, the bandwidth was approximately 
10% of the observing frequency, and the resulting 1% second- 
order correction would be of order a radian for typical val- 
ues of total electron content for the ionosphere. We fit the 
amplitude-weighted phase averages calculated for each chan- 
nel to a quadratic polynomial and used this polynomial to cal- 
culate the appropriate phase to subtract from the visibilities in 
each channel. This yielded the amplitude as a function of fre- 
quency, as shown for example in Figure [2 for the short Fort 
Davis-Pie Town baseline. 

To verify that our calibration procedure did not introduce 
any unintended biases, we calculated two-dimensional autocor- 
relation functions (ACF) for each data set using the visibilities 
from a calibrator source. As expected, we obtained results flat 
in both time and frequency lags, except for a spike at small time 
and frequency lags due to the effects of noise, broadened from 
a delta-function by smoothing. 

2.2. Observational Results 
2.2.1. Dynamic Spectra 

Initially, we had hoped to measure the angular broadening of 
the pulsar by measuring the phase variations of the diffrac- 
tion pattern; however, several complications derailed our ef- 
forts. Most significantly, the scintillation bandwidth was much 
larger than expected. As the phase variations of the diffraction 
pattern have this characteristic bandwidth (Desai et al. [1992 1, 



they could not be distinguished from ionospheric phase, phase 
slope, and curvature. Moreover, the angular scale inferred from 
the scintillation bandwidth is so small that we would not expect 
to detect phase variations on any Earth-based baseline, at 327 
MHz or higher observing frequency. 

We chose to focus on the Fort Davis-Pie Town baseline, 
which had the least interference and best signal-to-noise ratio. 
Figure [2 shows plots of the visibility amplitudes for the pul- 
sar on this baseline for both epochs. The white vertical lines 
correspond to the channels omitted at the ends of each IF, and 
the white horizontal lines correspond to breaks in the scan or 
to time ranges flagged because of RFI. At both epochs, we 
detected scintillation maxima that typically persist for ten to 
twenty minutes and that span several IFs, suggesting a charac- 
teristic width in frequency on the order of ten megahertz. 

2.2.2. Scintillation Bandwidth and Timescale 

We formed autocorrelation functions of the data shown in 
Figure ^ Figure |2 shows an example, a composite autocor- 
relation function formed for all the data. We then fit models 
to the autocorrelation functions, using the form expected for 
Gaussian spectra: a Lorentzian for differences in frequency and 
a Gaussian for differences in time (see Gwinn et al. 1998 1. 
We found values for the decorrelation bandwidth of Av = 
14.7 + 0.5 MHz for the first epoch and Av = 16.5 + 0.5 MHz 
for the second, measured as the half-width at half maximum 
of the autocorrelation function in frequency. A fit to a com- 
posite autocorrelation function, using all of our data, yielded 
Av = 15.7 + 0.2 MHz. This composite function is the sum 
of the autocorrelation functions for the two epochs, weighted 
by the number of samples. The quoted errors are standard er- 
rors for the fits; we discuss the expected uncertainties further 
below. The scintillation timescale tjss was about 1000 sec, as 
discussed below. TableHsummarizes the results. 

As Figure ^ shows, the number of independent samples 
with dimensions Av x tiss is small, so that the accuracy of 
our measurement is limited by number of samples. Phillips & 
Clegg ( 1992) argue that the fractional accuracy of such a mea- 
surement is given by the number of independent samples, so 
that 

cr(Av) 2 /Av = l/(W,xJV/) = tiss/T x Av/B. (1) 

For our observations, N, = T/t;ss ~ 14, and Nj = B/Av ss 2. 
We thus expect cr(Av) * (1/ V28)Av ~ 3 MHz. 

Because the total bandwidth is only twice the decorrela- 
tion bandwidth, we used Monte Carlo simulations with a range 
of parameters to estimate errors as well. Using the stationary- 
phase technique (Gwinn et al. U998> . we formed 1500 dynamic 
spectra with the same N t and Nf as our observations, with 
"true" decorrelation bandwidth evenly distributed between 
and 50 MHz. We then found the autocorrelation function for 
these synthetic spectra, and fitted for the decorrelation band- 
width, as described above for the actual data. We then found the 
distribution of "true" decorrelation bandwidth that contributed 
to Monte Carlo results close to our measured value. Thus, we 
inferred a parent distribution for our observed value. 
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We found that the parent distribution was sharply peaked 
near the measured value of 16 MHz, but had an extended tail 
toward larger values. Of those spectra that yielded values near 
our measured value, 67% came from the interval between 13 
MHz and 23 MHz. The endpoints of this interval have about 
equal probability, according to our Monte Carlo simulations; 
their asymmetry about the most probable value of 16 MHz re- 
flects the tail toward larger Av. The fact that we obtain values 
close to 16 MHz independently, from both of our observing 
epochs, provides additional evidence that 16 MHz is the "true" 
underlying value; we did not demand that our Monte Carlo sim- 
ulations reproduce this additional fact, to keep them manage- 
ably simple. This simulation suggests that the uncertainty sug- 
gested by Eq.^is approximately correct, with additional prob- 
ability on the side toward larger Av. The last line in Table \I\ 
indicates our adopted values for uncertainties, from the Monte 
Carlo simulations, of 16^ MHz. 

The Monte Carlo simulations indicate that the actual decor- 
relation bandwidth is unlikely to be about 1 MHz, as would be 
expected from some previous measurements (Paper 2). To test 
this hypothesis, we formed 1000 spectra with true decorrela- 
tion bandwidth of 1 MHz, and then found the autocorrelation 
function and fitted value for the decorrelation bandwidth. None 
approached 16 MHz. The largest value obtained was 4.7 MHz. 
We conclude that the probability of the decorrelation band- 
width actually being close to 1 MHz is much less than 0.1%. 

Similarly, we measured the scintillation timescale as the 
half-width at half-maximum of the time correlation function 
with time. The measured scintillation time scales were tiss - 
885 ± 9 s for the first epoch and t !SS = 1630 + 80 s for the 
second. The composite autocorrelation function yielded a time 
scale of tjss = 1020 + 70 s. The expression of Phillips & Clegg 
( 1992 1 would indicate that the true value was 1000 + 200 s. 
Monte Carlo simulations of the parent distribution indicate a 
confidence interval of tiss = 1000*™ s. The larger estimated 
errors we find from Monte Carlo simulations, for both tiss an d 
Av, may result from departures of the correlation function from 
its expected form, for small N, and Nf. We adopt the central 
value from the composite autocorrelation function in time, and 
the uncertainty from Monte Carlo simulations or the expression 
of Phillips & Clegg, as shown by the last line in Tabled 

2.2.3. Fine-Scale Scintillation 

A finer-scale variation within the scintillation maxima is appar- 
ent within the large-scale structures in Figure^ This variation 
is not noise: it is a real modulation of intensity. Figure[3]shows 
an expanded view of part of a scintillation maximum in Figure 
^ These finer-scale scintillations are also visible in the autocor- 
relation function, when the domain of integration is limited to 
a single peak of the scintillation spectrum in Figure^ Analysis 
indicates that they have bandwidth Av w 0.5 MHz and time 
scales of 1 to 3 minutes. The modulation index of this finer- 
scale structure is ntf » 0.2. These values are in approximate 
agreement with the scintillation parameters other groups have 
reported for PSR J0437-4715, as we discuss further in Paper 
2. 



For the fine-scale scintillations, the autocorrelation func- 
tion shows differences from the forms described in 32.2.21 
above, even when the reduced modulation index is taken into 
account. These might arise from the modulation introduced by 
the larger-scale scintillation; or by variations in the characteris- 
tic scales of the fine-scale scintillation, perhaps similar to those 
described by by Gothoskar & Gupta ( 2000 ), but within a single 
observation. The number of independent samples of the fine- 
scale scintillation is large enough to have negligible contribu- 
tion to the expected errors. We adopt the standard errors for our 
fits, as our best estimates. Variations of the properties of small- 
scale scintillation with epoch, as Gothoskar & Gupta observe, 
could be greater than these errors, but we cannot assess them 
from our data. In practice, longer and greater-bandwidth ob- 
servations of this pulsar, or observations at lower frequencies 
where scales are smaller, are probably required to understand 
statistics of the finer-scale scintillation. We discuss a possible 
explanation of the fine- and wide-scale scintillation in Paper 2. 

3. Discussion 

3.1. Scintillation Velocity 

The combination of scintillation bandwidth and time scale 
yields an estimated velocity for the scintillation pattern, usually 
dominated by the velocity of the pulsar. The relation (Gupta et 
al. lT994l : 

VAv/ss DX 

V ISS =A V -L-^ (2) 

/ tiss 

connects scintillation velocity Viss to the observables Aviss 
and tiss, along with the distance of the pulsar D, a constant 
Ay, observing frequency /, and X, the ratio of the distance from 
observer to scattering material to that from scattering material 
to pulsar. For material halfway X = 1 ; for material close to the 
observer X < 1 . 

The calculated scintillation velocity Vjss varies along 
with the scintillation bandwidth. Gothoskar et al. (|2000l and 
Nicastro et al. ([1995 1 both report variation of the inferred Viss 
by about a factor of 2, as inferred from individual observa- 
tions. The Viss obtained for our adopted values for the wide- 
band scintillation in Tabled with A v - 3.85 x 10 4 km s~ ! x 
(1 ■ GHz 1 sec)(l MHz • lkpc) 1/2 (Gupta et al. 113531 . and an 
assumed distance of 150 pc (van Straten 2001 1, and X — 
1, is 180^y kms~'. This is close to the average value of 
170 km s~ reported by Johnston et al. ( 1998 1, and consistent 
with the 231 km s _1 reported by Gothoskar & Gupta (2000 1. 
The proper-motion velocity is 100 km s (van Straten et al. 
2001 1. The values calculated from our two epochs separately, 
200 km s _1 and 1 10 km s"\ differ with marginal significance 
from the adopted value; more extensive observations could 
reveal whether the time-variability observed by Gothoskar 
& Gupta holds for the wide-band scintillation. Gothoskar & 
Gupta ( 2000 1 inferred from the relatively high scintillation ve- 
locity that the fractional distance of the scattering material, x 
(see Eq.|2), was larger than 1, and surmised that the scattering 
material likely lies at the surface of the Local Bubble. As dis- 
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cussed in Paper 2, some evidence suggests that the scattering 
lies in a thin layer, w 10 pc from the Sun. 

3.2. Comparison with Other Pulsars 

PSR J0437-4715 has among the lowest dispersion measures 
of strong pulsars, and it is perhaps unsurprising that the scin- 
tillation bandwidth is among the broadest observed, scaled to 
a single frequency. Comparisons with other pulsars are usually 
done by the typical pulse broadening time r (actually the time 
required for a broadened pulse to reqch 1/e of its maximum) 
(Pynzar' & Shishov ll997> . The scintillation bandwidth is re- 
lated to the typical pulse broadening time, r, by the equation: 

Av ISS T = C 1 /(2n) (3) 

where C\ is a constant of order unity (see Taylor et al. 1993 
Lambert & Rickett lT999t . We follow Pynzar' & Shishov ( fT997b 
and adopt C\ = 1.0, and scale to observing frequency /o = 330 
MHz using r oc /~ 4 - 4 . Figure 0] compares the scintillation pa- 
rameters of PSR J0437-4715 with those of other pulsars. The 
broader-Av/ SS scattering of PSR J0437-4715 falls below the 
scaling found by Pynzar' & Shishov (1997) for nearby pul- 
sars, as does the scattering of PSR B0950+08 (Phillips & Clegg 
I1992> . Interestingly, some earlier measurements of the scintil- 
lation properties of PSR B0950+08 found spectral structures 
that would indicate much larger r, placing it near the line ex- 
trapolated from other nearby pulsars. Those measurements also 
indicated frequency and time structures unexpected for scin- 
tillation, and large variations with observing epoch (Roberts 
& Abies 1982 1. The fine-scale structure we observe for PSR 
J0437-4175 lies above this extrapolated line, although obser- 
vations of this object by others yield values close to the line. 
We discuss observations of PSR J0437-4715 by other groups, 
and present a possible interpretation, in Paper 2. 
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4. Conclusions 

We find two scales of scattering for the nearby pulsar PSR 
J0437-4715. The broader scale has scintillation bandwidth 
Av ISS = 16+^ MHz, and time scale t !SS = IOOO^q sec. The 
narrower scale has Avjss = 0.5 + 0.1 MHz, and time scale 
tiss - 90 ± 20 sec. Both yield scintillation velocities close to 
the proper-motion velocity found by van Straten (2001 1. The 
pulse broadening time inferred from the broader scintillation 
bandwidth is less than that expected on the basis of observa- 
tions of other nearby pulsars. 
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Fig. 1. Visibility amplitude as a function of frequency and 
time for observations of PSR J0437-4715 on the Fort Davis- 
Pie Town baseline on 5 November 1996 (bottom) and 8 April 
1999 (top). 
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Fig. 2. The composite ACF obtained by suitably combining the 
ACFs from the two observations. The best-fitting Lorentzian 
has a characteristic width of Avjss = 15.7 ± 0.2 MHz. 
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Fig. 3. Variations in the visibility amplitude within a scintilla- 
tion maximum appear to occur over a frequency scale on the 
order of 1 MHz. 
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Fig. 4. Pulsar pulse broadening time plotted with dispersion 
measure, adapted from Pynzar' & Shishov ( 1997 1. As they 
discuss, the solid line shows the scaling of pulse broaden- 
ing with distance expected for a uniform medium. Stars show 
the wide-band (4-pointed) and fine-scale (5-pointed) values we 
measured for PSR J0437-4715; error bars and variation be- 
tween epochs lie within the symbols. Triangle shows value for 
PSR B0950+08 measured by Philli ps & Cl egg (TT9921 Data for 
other pulsars from Cordes & Lazio (2002i. 
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Table 1. Measured Scintillation Parameters for PSR 
J0437-4715 





Wide-Bandwidth 


Narrow-B andwidth 


Epoch 




tiss 




tiss 


(MJD) 


(MHz) 


(sec) 


(MHz) 


(sec) 


50392 


14.7 


885 


0.6 


90 


51276 


16.5 


1630 


0.4 


90 


Composite 


15.7 


1020 


0.5 


90 


Adopted 


16+* " 


1000+^ 8 


0.5 ±0.1 " 


90 ± 20 * 



8 Error estimated from Monte Carlo simulation. * Error estimated from 
standard error from fits; see 8 12.2.31 



